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Abstract 

During catastrophic processes of environmental variations of a thermody- 
namic system, such as rapid temperature decreasing, many novel and complex 
patterns often form. To understand such phenomena, a general mechanism 
is proposed based on the competition between heat transfer and conversion 
of heat to other energy forms. We apply it to the smectic-A filament growth 
process during quench-induced isotropic to smectic-A phase transition. Ana- 
lytical forms for the buckling patterns are derived and we find good agreement 
with experimental observation [Phys. Rev. E55 (1997) 1655]. The present 
work strongly indicates that rapid cooling will lead to structural transitions 
in the smectic-A filament at the molecular level to optimize heat conversion. 
The force associated with this pattern formation process is estimated to be 
in the order of 10 -1 piconewton. 
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When a thermodynamic system in some initial stable phase is subjected to rapid or 
abrupt environmental variations such as decreased temperature, varied ionic conditions, 
etc., many novel patterns will often form. For example, in volcanic eruptions streams of lava 
are spurted out and cooled rapidly, during which many tensile bubbles are spontaneously 
formed and the lava becomes porous material [1, 2]. Similar bubbles and strips are seen 
in supercooled biomembranes [3]. And in rapid growth process of carbon nanotubes initial 
straight tubes are observed to change into regular helical shapes, even at fixed environmental 
conditions [4]. Other puzzling phenomena include the kink instabilities of short DNA rings 
in solution: on adding some specific ions, a circular ring changes abruptly into polygonal 
ones [5]. 

All these transitions and pattern formations are radically different from those slow, quasi- 
static processes discussed in conventional textbooks, such as the isothermal solid-liquid 
transitions. In the solidification process of water, for example, the transition takes place 
very slowly and at any time during this process, solid and liquid phases can be regarded 
as in equilibrium with each other and with the environment, the bulk temperature is kept 
constant and latent heat of transition is released. However, for those catastrophic processes 
mentioned above [1-5], the quasi-static and quasi-equilibrium assumptions are no longer 
satisfied since the system is now far away from equilibrium and changes rapidly over time 
[6]. Therefore, we have to create new ways and new insights to properly understand them. 
A deep understanding of such phenomena is needed both for its own sake and for various 
practical purposes, such as the optimal growth conditions for some important materials. 
In recent years many research efforts have been devoted to this field [6-10]. For example, 
in Refs. [6] and [7] (and in other references) the spinodal decomposition and formation 
of separate patterns in a binary mixture is studied both theoretically and experimentally, 
using respectively linear stability analysis and x-ray scattering technique. And Refs. [8-10] 
and [6] have investigated the solid-liquid interphase patterns during directional solidification 
processes induced by large temperature gradient. A very important aspect of the mechanism 
for these pattern formation phenomena is that the pattern formation is controlled mainly by 
the relatively slow diffusive transport process of heat energy out of the system [6]. Thus, It 
seems that this mechanism can not give a satisfactory explanation to those pattern formation 
processes mentioned in the first paragraph of this paper. This is partly because of the fact 
that in these systems the effects of other energy forms become comparable to that of heat 
energy [1-5], and they are anticipated to take a significant role in determining the formed 
patterns. 

In this paper, we try to understand the nature of these rapid pattern formation processes 
by regarding them as quench-like ones induced by, say, temperature decreasing. Our main 
reasoning is as follows. In a quench process caused by rapid cooling, the thermal energy or 
heat of the system is required to be evacuated as quickly as possible. However, this demand 
can not be met just only through the relatively slow heat exchange process occurred at the 
system surface, and a great part of the thermal energy are thus anticipated to be converted 
into other energy forms [11, 12]: mechanical, elastical, and maybe even electromagnetical 
ones, which can accommodate large amount of energy quickly and efficiently Therefore, 
during the quench process, because of the competition between heat transfer and heat con- 
version, the structures of the system, both at macroscopic level and at molecular level, might 
change dramatically, and many novel and complex patterns are anticipated to form sponta- 
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neously. It resembles in some sense the rapid expansion of ideal gases, which is regarded as 
an adiabatic process since the system has no time to exchange heat with the environment. 

To check the validity of this perspective and show how to implement it in actual practice, 
we investigate into a specific phenomenon based on this insight, namely, the growth and 
pattern formation process of smectic-A (Sm-A) filaments from an isotropic (I) liquid phase 
induced by rapid cooling [13]. Although the subject is quite specific, the general reasoning 
is the same and can be applied to many other pattern formation processes. Our theoretical 
results show that: On cooling, an initial straight filament will buckle and take on some 
specific kind of curved configurations which can be well described by elliptic integrals. Good 
agreements with the experiment of Ref. [14] have been attained, indicating the usefulness and 
robustness of our treatment. As also speculated in Ref. [14], the present work confirms that 
structural transitions in the Sm-A filament at the molecular level are required in order to 
optimize heat conversion. The force associated with this pattern formation is also discussed 
in this paper. 

As observed in the I-Sm-A phase transition process [14] , at the initial stage, thin, hollow, 
and straight Sm-A liquid-crystalline tubes appear from the I phase on cooling. As the 
temperature is further decreased, these straight filaments buckle at the growing ends and 
take on somewhat complex serpentine forms, the amplitudes of which become more and more 
large. The filaments are metastable and eventually transform to compact domains after the 
filaments become extremely convoluted. Although this transition phenomenon seems to be 
very complex, we found its physics is actually simple and the formed patterns can be well 
described by analytical functions. 

As we have known, the net difference in the energy between a thin Sm-A filament and 
the I phase is the sum of the following three terms [13, 14]: 

(1) the volume free energy change due to I-Sm-A transition F v = —g V = —n{p 2 — 
Pi) 1 9ods, where go is the difference in the Gibbs free energy densities between I and Sm-A 
phases and V is the volume of the Sm-A nucleus, and p Q and pi are respectively the outer 
and inner radii of the Sm-A microtube, s is the arc length along the central line r of it; 

(2) the surface energy of the outer and inner Sm-A-I interfaces Fa = 7(A> + AA = 
2vr(p + Pi) 7 / ds, where 7 is the Sm-A-I interfacial tension, and A Q and Aj, are the surface 
areas of the outer and inner surfaces, respectively; 

(3) the curvature elastic energy of the Sm-A filament F c = (ku/2) /(V • N) 2 dV + (2kn + 
kz>)(p — Pi) § KdAi — k c J n(s) 2 ds + nkn ln(p /pj) / ds, here k^ is the Oseen- Frank elastic 
constants and k 5 = 2k 13 — k 22 — k 2i , k c = nku{p 2 — pf)/4, and N and K are respectively 
the director and Gaussian curvature defined on the inner surface of the Sm-A tube; k is the 
curvature defined along r [14, 15]. 

The total energy difference is thus F = Fy + Fa + F c = k c J k,(s) 2 + /(A + F\)ds, where 
A = -nku\xi(p / p-i) + 2n^(p + p.j) and n = —ngo(p 2 —p 2 ) [14]. We call F the shape formation 
energy of the Sm-A filament. The variation equation 5F = yields the equilibrium-shape 
equation of the filament: 

k c {^ - 2kt 2 + 2k) - (A + n)/c = 0, (1) 
Akr + 2kt = 0, (2) 

where r is the torsion of the filament's central line r, according to the Frenet formula [16]; ( ) 
and ( ) mean respectively first and second order differentiations with respect to arc length 
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s. 

It is obvious that a straight line is always a solution of the above equation, since its 
k = t = 0. The corresponding shape formation energy of a straight filament is F = (A + n)Z, 
where / is the length of the straight filament. In the rapid cooling process, it is impossible 
to effectively transfer redundant heat out of the system by means of heat exchange with 
environment, and Sm-A straight filament will form spontaneously and quickly to partly 
accommodate the redundant energy, as long as A + n > 0. Straight Sm-A filaments have 
been observed to form at the initial stage of the growth process [14]. However, since the 
value of g increases with temperature decreasing [14, 17], while 7 exhibits no temperature 
dependence or a weak one [18], A+n will decrease with cooling, and the equilibrium threshold 
condition of F = yields the criteria for the growth of a straight filament as 

nkn \n(p /pi) + 2n>y(p + p t ) - ng (p 2 o - p 2 ) = 0, (3) 

beyond which, forming a straight filament will no longer help to absorb thermal energy. This 
situation is unfavorable to the quench process, and consequently the resultant remnant part 
of energy will prevent the straight filament from keeping stable. Then a shape deformation 
will be induced, which would lead to another solution of the shape equation to convert as 
much amount of thermal energy into elastic deformation energy as possible. In other words, 
the straight filament will buckle to a curved configuration which, on the one hand, its k and 
r should satisfy (1) and (2) so the new pattern will be stable, and on the other hand, the 
stable new configuration should be capable of absorbing as much thermal energy as possible. 
Then, to what configurations will the filament deform? 

Now we would prove that the observed Sm-A serpentine planar patterns shown in Ref. [14] 
are just the allowed solutions of Eqs. (1) and (2). For the convenience of comparison with 
experimentally observed planar curved patterns, we will confine ourselves to the planar case, 
where r = 0. Then Eq. (2) is automatically satisfied and (1) reduces to 

k c (n 4 + 2nk) - (A + H)k 2 = 0, (4) 

here, we have multiplied both sides of Eq. (1) with k. Define the filament central line tangent 
unit vector to be t(s) = (cos0(s), sin0(s), 0), then k 2 = <f) 2 . Inserting this into (4), we know 
that is determined by the following equation: 

4> 2 = -2t7 + /-cos(0-0 o ), (5) 

where 77 = —(A + n)/2k c [14] and / is a positive constant, (later, we will show that / 
is actually the mean repulsive force exerted on the growing end of the filament), 0o is 
an integration constant. The solution of (5) is a periodic function, its form in the range 
(j) G (0, arccos(l — 2m)) is expressed as [14] 

\2k c . /l — cos0 , — . , . 

y— F(arcsmy — — — , y/m) = {s - s ), (6) 

here, 

m=l-^€(0, 1) (7) 
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and So is another integration constant, F(ip, k) is the first kind elliptic integral. In the same 
range the filament shape is determined, up to an additive constant vector, by [19] 




1 — cos0 — /l — cos0 — / ; — 

z£/(arcsm \ , y/m) — Ffarcsm \ , v m J? — \/2m + cos0 — 

V 2m V 2m v 

(8) 

where E(ip, k) is the second kind elliptic integral. The trajectory of the filament in the whole 
range can also be obtained easily from (5), (6) and (8). The general shapes are shown in 
Fig. 1 in the cases ofm<l/2or^>0 (dotted curve), m = l/2or^ = (dashed curve) 
and m > 1/2 or rj < (solid curve) for one period of the trajectory. The total arc length of 
each period is calculated to be [19] 



P = ^K{y/m) (9) 



and the shape formation energy density of each period is [19] 

* = 4 - 2r >) ds = f [ 2L (vM - 3 + 4ml , (10) 

i Jo 1 J 

where K{k) and E(k) are respectively the first and second kind complete elliptic integrals, 
and L(v/m) = E(y/m)/ 'K(y/m). 

At the initial stage of buckling, the filament will deform to a configuration to make 
the threshold condition x = be satisfied, which leads to m = m^ = 0.339175 at this 
stage. The corresponding deformed shape is the dotted curve of Fig. 1, it is very similar 
to the observed initial buckling shape Fig. lb and lc of Ref. [14]. Furthermore, for this 



shape, the amplitude of the deformation is calculated to be Y rnax = yAmk c /f from (8), 
and the distance traversed along the growth direction for one half of the trajectory period is 



Xmax = y 2k c / f (2E(^/m) + K (^/m) , thus Y max /X max = 0.737 for the above mentioned value 
of m,j. On the experimental side, we can estimate the values of this ratio to be respectively 
about 0.85 and 0.64 from Fig.lc and Id of Ref. [14], and the mean value is 0.745, all in close 
agreement with the theoretical value. 

As mentioned before, after buckling, the filament will try its best in adjusting its struc- 
tures to accommodate thermal energy. Differentiation of x lYy (10) with respect to m, we 
find that the maximum value of the shape formation energy density x is reached at the point 
where m = 0.9689, and in the range between m,j and m Q , x is a continuously increasing 
function of m. Thus, after the buckling process has been triggered, m will try to increase 
from its initial value m«, till the value m is reached. Since m is determined by (7), we antic- 
ipated that A + n will change from negative value ( and hence m < 1/2) to zero (m = 1/2) 
and then to become more and more positive (m > 1/2). The question is: How can this be 
achieved by the Sm-A filament? 

We know that go increases with temperature decreasing, and 7 exhibits little temperature 
dependence, and that during the cooling process there is no significant changes in the filament 
radii p Q and p^, therefore, the only possible reason accounting for A+n changing from negative 
to positive values is that k u , the quantity related to the molecular structures of Sm-A 
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filaments, rapidly increases its value. In other words, quench may induce some structural 
transition in the Sm-A filament at the molecular level. (In Ref. [14] it is speculated that 
such a transition is the trans to cis transition occurred in the alkyl chains of the Sm-A 
liquid crystal material. Further experimental observations are still needed in this respect). 
Consequently, the initial curved pattern will take the form demonstrated by the dotted curve 
in Fig. 1, later formed segments will take the form of the dashed curve, and still later formed 
segments will be the form of the solid curve. This trend is demonstrated more clearly in 
Fig. 2. Buckled pattern resembling that of Fig. 2 has indeed been observed and recorded 
in the experiment of Ref. [14] (see Fig. 1 of [14], especially Id), showing the correctness 
of the pattern formation mechanism based on competition between heat transfer and heat 
conversion. 

When m is increased to the value m c = 0.731183 < m , different segments of the same 
filament begin to intersect, as shown in Fig. 3. At this stage, the formed serpentine pat- 
terns may begin to collapse because of possible strong contacting forces, and eventually the 
filament will transform to compact domains [14]. Thus actually the optimal value m can 
not be reached by the system. 

So far, we have not discussed the physical meaning of the parameter / in Eq. (5). Now 
we begin to investigate into this question. Our result is that / is actually the mean repulsive 
force exerted on the filament growing end. There has been clear experimental evidence that 
both the assembly and disassembly of microtubular filaments can generate force; however, 
only limited quantitative data are available on the actual magnitude of these forces [20] . In 
the case of Sm-A filaments discussed here, we estimate the force exerted on the ends to be the 
order of 10 _1 piconewton. This magnitude is reasonable as indicated by various experiments 
on thin filaments [20, 21]. The forces associated with the assembly and disassembly of 
filaments may be of vital biological significance. For example, it has long been speculated 
that the motion of chromosomes during mitosis of the cell cycle is caused by the assembly 
and disassembly of cytoskeletal microtubes [20]. 

To show that the parameter / in Eq. (5) corresponds to the average force exerted at 
the ends, we will rederive (5), the stable planar shape equation, from another way. The 
already grown Sm-A filament can be viewed as an inextensible string. When subjected 
with the constraint of fixed end-to-end distance, the total shape energy is expressed as 
F = J k c t 2 + A + n — f • t ds, with t(s) = (sin 9 cos 0, sin 9 sin 0, cos 0) being the string 
direction vector; here f is a Lagrange multiplier corresponding to this constraint, which, 
according to thermodynamic principles, is just the average force exerted by the environment 
on the filament ends. The first variation 5F = leads to the shape equation expressed in 
terms of 9, <fi and their differentiations with respect to s, For our planar case with 9 = %/2, 
the non-straight stable configurations are found to be determined by the following equation 

2 = C + fflcos(0-0[, + 7r), (11) 

K c 

where <p' is the angle between the x-axis and the direction of the average force f. From 
(11) we know the filament tangent direction (growth direction) t is on average antiparallel 
to f, therefore f must be a repulsive force. Comparing Eq. (11) with (5), we can obviously 
infer that the parameter / in (5) is just the magnitude of the average force f and C = —2r/ 
in (11). This force is induced by thermal current, and we have assumed it to be constant 
throughout the growing process. 
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The magnitude of the force can also be estimated from the buckling amplitude of the 
filament. As mentioned before, the deformation amplitude Y max = ^Amk c / f, therefore 
/ = Amkc/Y^^. At the initial stage of buckling, m = = 0.339175, we take ku = 10~ 6 
dyn according to Ref. [13], and estimate that p Q = 1.3 pm, = 0.8 /iin, Y max = 10.4 pm 
from fig.lc of Ref. [14]. Thus, the average repulsive force is calculated to be / ~ 1.4 x 10~ 13 N 
or 0.14 pN. It is interesting to note that the forces related to polymeric strings are often in 
the order of 10~ 2 to 10 1 pN [20, 21]. 

In summary, in this paper we have proposed a theory to explain the novel pattern forma- 
tion phenomena induced by catastrophic variations of environmental temperature or other 
parameters. We then applied this general insight in understanding the pattern formation 
process of smectic-A filaments grown from an isotropic liquid phase. Exact buckling patterns 
have been obtained and compared with experiment, and we achieved excellent agreements 
with experiment. We also suggested that temperature decreasing leads to some structural 
transitions at the molecular level for the Sm-A filament, possibly the trans to cis transition 
of the alkyl chains of the material used in the experiment of [14]. The force associated with 
the pattern formation process of smectic-A filaments is estimated to be in the order of 0.1 
pN. 

One of the author (Z.H.) appreciates valuable discussions with Drs. Yan Jie, Liu Quanhui 
and Zhao Wei. We are also indebted to the anonymous Referee, who makes us take notice of 
Refs. [6] and [7] . This work is partly supported by the National Natural Science Foundation 
of China. 
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FIGURES 



FIG. 1. Shapes of one period of the buckled patterns (Eq. (8)) for m = 0.339175 (dotted line), 
m = 0.5 (dashed line) and m = 0.7 (solid line). The length unit is set to be y/2k c / f. 

FIG. 2. Simulated buckling filament shape during the growth process. Later formed segments 
have longer period arc length and larger amplitude, caused by increased values of k\\. The unit 
length is set to be J2k l c / f, where k l c is the value of k c at the onset of buckling. 



FIG. 3. Onset of self-intersection of the buckled pattern at m = 0.731183. The appearance of 
self- intersections might trigger the transition from filament structures to compact domains. The 
length unit is the same as in Fig. 1. 
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